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1. Introduction 

The problem of bound states in the relativistic quantum field theory with 
large coupling constant is an extremely difficult problem. To understand 
the structure of the bound state is a very ambitious goal which assumes 
the solution of a whole spectrum of tightly connected problems, such as 
confinement, chiral symmetry breaking phenomenon, and many, many others 
which are greatly important in the low energy region. 

A less ambitious purpose is the study of the hadron wave function (wf) 
with a minimal number of constituents. As is known such a function gives 
parametrically leading contributions to hard exclusive processes. In this case 
the quark and antiquark are produced at small distances z ~ 1/Q -^ 0, 
where Q is the typical large momentum transfer. Thus, we can neglect the 
z'^ dependence of the wave function of the meson with momentum p and can 
concentrate on the variable zp ~ 1 which is order of one. Therefore, the 
problem is drastically simplified in the asymptotic limit and we end up with 
the light cone wave function, (f){zp, z'^ = 0). 

The corresponding wave functions have been introduced to the theory 
in the late seventies and early eighties [0] in order to describe the exclusive 
processes in QCD. We refer to the review papers 0,0, on this subject for 
the detail definitions and discussions in the given context, but here we want 
to make a short remark concerning the very unusual properties of the non- 
perturbative light-cone wave functions of the leading twist. The information 
which can be extracted from the QCD sum rules method 0,0 unambigu- 
ously shows the asymmetric form of the distribution amplitudes and this 
property was unexpected and even suspicious to many physicists. 

However, to study the fine aspects of the theory (see explanations below) 
of the exclusive processes and to extend the area of applications, we need 
to know not only the dependence wf on the longitudinal variable Xi, but on 
transverse variable k'j_ as well. In particular, as is known |0],0, the Sudakov 
suppression should be taken into account in order to integrate correctly over 
the " endpoint region" [xi -^ 1). The dependence of wf on the transverse 
momentum /c^ plays an important role in such calculation. Besides that, the 
transverse size dependence plays a key role in color transparency physics. 

• The main goal of the paper is the calculation of the few lowest moments 
(fc^") of the transverse distribution. Besides that we formulate the definition 
of nonperturbative wave function and find some constraints on it. In partic- 
ular, we analyze the asymptotically distant terms {k"^), n — > oo in order 
to model the ip{k\,Xi). The byproduct of this consideration is analysis of 
the factorizability of ilj{k\,Xi) = cl){xi)ip{k'j_), which it turns out, does not 
work. Finally, we model wf which meets all constraints. Analogous problem 
for the longitudinal distribution (^") has been formulated and solved within 
QCD sum rules method in ref. 0. 

The paper is organized as follows. In the next section we derive the 
explicit relations between two, apparently different, values: 
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l.{k'^)p for pseudoscalar wf, (pp = {0\d{z)i^5u{—z)\7i) on the one side and 
2. vacuum expectation values (VEVs) (qC^^q) on the other side. Such 
exact (in chiral hniit) connection is not very surprising because vr meson is 
the Goldstone particle, strongly interacting with vacuum fluctuations. Well- 
known example of such connection is, of course, the formula 

(0|Ji75u|7r) = —-{0\dd + uu\0), (1) 

which relates vr meson matrix element and some vacuum characteristic. 

In the section 3 we derive the analogous relations for {k'^)A for axial 
(leading twist) wf. In this case the corresponding formulae are not exact 
even in the chiral limit, but we argue that their accuracy is very high and 
corrections are of order as/vr ^ 10%. We will see that for both cases (axial 
and pseudoscalar ) the problem is reduced to the evaluation of the mixed 
vacuum condensates {qG^^q). 

In order to estimate them we consider in section 4 some special kind of 
sum rules which are very sensitive to the VEVs we are interested in. We find 
a large violation ( order of factor 3) of factorization for mixed condensates 
{'ii.^ ij.uG pLvY l) of dimension seven. It automatically leads to a noticeable 
fluctuations of the momentum {k'\)A > (^Di and leads to the well-spread 
out wf (in the momentum space). Let us note, that if the factorization 
would work, we would get {k\)A,p ~ {k±)'pA- Such a relation gives the 
qualitatively different behavior for wf and means a strong concentration of 
the distribution density around point k'j_ = 0. 

In the section we derive an "almost" exact relation between mixed vacuum 
condensates {q{igo^\aG\„Yq) = {rnQ)"-{qq) for arbitrary n. 

The section 5 is the main part of this paper. We formulate the definition 
of nonperturbative wave function and give some constraints on it. Finally 
we model it. 

2. The QCD vacuum condensates and {k'^)p- 

We define the pion pseudoscalar wave function {0\d{z)i-y5u{— z)\n{q)) = 
~ f'^' (f)p and the corresponding mean values of the quark transverse distri- 
bution in the following way: 

{0\d^^,i^D,t,ruHq)) = -^(-^')"^^J^(^1")p> (2) 

where iD^ = id^ + gA"^^ is the covariant derivative and a transverse vector 
t^ is perpendicular to the hadron momentum q^. The factor ^ ^2ny\" ^^ related 
to the integration over angle in the transverse plane: / t/0(cos 0)^"/ J d(j) = 
(2n-l)!!/(2n)!!. 

We removed the common factor — ^ in order to reproduce the matrix 

element (|I|) without derivatives. To find {k'j_)p we follow the paper [10| and 



consider the following matrix element 

{Q\di^^iD^iD,u\Tx{q)) = ^j^Ag^,, {kl)p = -2A, (3) 

where we dropped off the kinematical structure g^gj, by the reason which will 
be clear soon. Multiplying (^) by g^^ and using the equations of motion and 
PCAC, we get (in the chiral limit, iriq ^ ): 

{0i^,iD,W,u\ix{q)) = -^AA = ^{qiga,,G%^q), (4) 

Jit J it ^ 

4(gg) 4 ^ 

where we use the standard value for parameter ml [^. It is clear that the 
skipped term is proportional to q^quQ^iv ~ "^^ ~ ttt^q and gives zero contribu- 
tion in the chiral limit. 

In order to find {k\)p we consider the following matrix element 



2{qq) 



{0\di'j5iD^^iDf,JD^JD^^u\7i{q)) = — (5) 



where we dropped off the terms which include the g^g^, structure by the 
reasons mentioned above. The PCAC leads to the following relation 

- ^ ^ ^ ^ 2 

{0\di-f5iD^,iD^,jD^.^iDf,^u\n{q)) = — ^(gP^iP/.2P/.3^M49), (6) 

Jtt 

where P^ = iD^ is hermitian operator and condensate {dP^j_-^P^^P^.j^P^^q) can 
be evaluated in a standard way. The result is: 

/r4\ -{q9^<yfiuGi,uCrxaGxaq) + 2{qg'^G^uG^yq) ,, 

**-^*" Trni ■ '" 

where G^^ = G^^^. 

The analogous relations for arbitrary n are drastically simplified in the 
limit when the space-time dimension d ^ oo. Of course, the real world 
corresponds to c? = 4, but we expect that the expansion 1/d^ works well 
enough and errors are of order 1/d ~ 25%. Anyhow, we are not going to use 
this limit for the numerical estimates. Instead, we want to use this limit in 
order to understand the general structure of the (k'j^) (behavior at ra — > cxd, in 
particular). In this limit, discarding the small terms related to the creation 
of the additional gg pair^j, we get the following simple formula for (fci")p 

^Indeed, the additional terms which appear in this derivation have the form 
[iD,a^j,vGfj,u\ ~ {iD^ji.G^j,v)^u ^ qju^'^qji'- It can be checked that these terms are neg- 
hgible smaU. 



expressed exclusively in terms of condensate {qiigcy^yGi^iuYl)^ which is the 
leading operator in d -^ oo limit: 

(«;i> =. c" <*2f^#f*::2> , d^oo (8) 

m) 

with some constant C. 

Now let us come back to the formula (^). It is clear that the calculation 
of the mean value of the quark transverse momentum in the pion and the 
evaluation of some vacuum characteristics are two sides of the same coin. 

Let us derive the same relations (^) in a different way. Consider for this 
purpose the following correlator 

T...^„ =1 j dxe'^%Q\T{di^,0^,..^Mx),u^,^,dm\^) = q^T^^n}iq') + -(9) 

with arbitrary local operator 0^j_^,,^^. We extracted the special kinematical 
structure q^T^^^y in such a way that the tenzor T{^^} does not depend on 
Qfj,, but depends only on metric g^^^y The leading contribution to T^n„} at 
q^ ^ oo can be easily calculated and is given by 

T{,M) = ^{QO,,..,„q) + O(^). (10) 

On the other hand, the only particle which can contribute to this correlator 
is TT meson, because a^ = u'julbd is conserved in the chiral limit and matrix 
elements (0|a^|7r') for the massive pseudoscalar particles are proportional to 
dutti, ~ rrig -^ 0. At the same time, the axial mesons do not contribute to 
the kinematical structure we chosen q^T^^^^. 

Thus, the whole result (^) comes exclusively from the vr meson and any 
corrections to this statement are proportional to the non-conservation of the 
current d^ai, ~ mq ; in particular, the next term proportional to 1/q^ comes 
with the small factor niq. Therefore, the correlation function T{^^} looks like 
a huge peak related to vr meson contribution plus some function on q^ which 
is nontrivial, but proportional to the small number mq. Using the dispersion 
relation and retaining only vr meson contribution to ImT^^^y we recover the 
formula (|^). The great usefulness of this derivation will be appreciated soon, 
in the analysis of the axial wf. 

3. {k^)A- Leading twist wave function. 

We define the pion axial wave function if-^q^cpA = (0|(i(-2)7^75u(— z)|7r(g)) 
and the corresponding mean values of the quark transverse distribution in 
the following way: 

{0\d^M^D,t,y-u\7r{q)) = ^/.g.(-t^)" ^ (2n)!! '" ^^^"^^^^ ^^^^ 



First of all let us show that the calculation of eq. (0) is reduced to the 
calculation of the matrix elements depending only on Gfiy. We start from 
the simplest case and consider the following matrix element p!0| : 

{^\hvl^l^[^D ^ll^D ^L2 + «^/.2^^m]^k(g)) = (12) 

From the definition of (A;^)^ it is clear, that {k'^)A = — 2A. At the same 
time, multiplying (|12D by metric tenzor and using equations of motion, the 
constant A can be expressed in terms of the following matrix element: 

{0\d'l,^l5i^-)(^\aGxaU\TTiq)) = —iUq„A (13) 

In order to express {k\)A_ in the same terms, we consider the following matrix 
element 

(0^7^75^ Y. («^m^^M2^^w^^/^4)wk(?))- (14) 

{|J'l,^J■2,tJ■3,^J.4} 

By the reason mentioned above, we have to consider on the right hand side of 
this equation the most general kinematical structure which includes only the 
first power of q^ and which is symmetric under /ii,/X2,/^3, ^^4 permutations; 
the rest indexes are cared by the metric tenzor g^^. Multiplying ( |14|) by g^i,, 



using the equation of motion and skipping all small terms which are related 
to creation of the additional qq pair (for justification, see below), we arrive 
to the following result: 

i(lpf^{k]_)A = g(0|rf7p75(^aAaG'Aa)^^|7r) + (15) 

3 - 1 - 

— (0|d7p755f^G'^^G'^^M|7r) + —(0 16/7^755-^ (G^pG^^ + G^,yG^p)u\7r). 

Again, the problem is reduced to the calculation of some matrix elements, 
which can be estimated by the method, described in the previous section and 
which we are going to use now. 

In order to calculate the {kj^A let us consider the following correlator 
(compare with (^): 



ij dxe"^^0\T{dj^j5tgaxaGxMx),^^75dm\0) = q^,T{q') + ... (16) 

In comparison with the previous case, the current J^ = djfj.'J^igcrxaGxaU is 
not conserved even in the chiral limit, but it is "almost" conserved (9^J^ ~ 
[iD^G^u) ~ qii^q in a sense, that the non-conservation is small and related 



to the production of the additional qq pair. From the point of view of the 
expansion T{q^) at large q^ it means that the main contribution is given by 

■7 

T{q^) = -i^iqigcrxaGxaq), (17) 

q^ 

like in the previous case ([TI1|); the corrections to this formula are suppressed 
by the loop factor — ^ 0.1 . We neglect them in the following. Using a 
dispersion relation and retaining the n meson contribution only Q, one gets 

(/cj^-gg ^_^^ - 3g -U.lGel/ . (18) 

The same procedure can be applied for any operator whose matrix element 
we are interested in. The result for {k'\_)A is: 

The problem, like in the previous case is reduced to the analysis of the mixed 
vacuum condensates. 

4. Analysis of the condensates. 

Let me, first of all, formulate the method we follow in this section in 
order to estimate the higher dimensional condensates. The idea is as follows. 
Consider the correlation function T^ = i J dxe*''^(T{ J^(x), J(0)}) at large g^. 
If the currents J^, J are chosen in such a way, that in the chiral limit the 
perturbative contribution is zero and the current J^ is "almost" conserved 
(in a sense of the previous section), we end up with the leading (at g^ -^ oo) 
contribution in the form ■^, like (0), plus some nontrivial function, but with 
small coefficient, like — , or rriq, in front of it. Besides that, if we know the tt- 
meson matrix elements (0| J^lvr), (0| J|7r) exactly, we can find the condensate 
(0) by collecting all leading terms proportional to 1/g^. To convince the 
reader in correctness of this method, we will explicitly calculate the loop 
corrections as well in the example which follows. 

In order to calculate the condensate {qg^<J^uGf^u<^xaGxo-q), let us consider 
the following correlator: 



tj dxe''^^{0\T{d^p^5^g(TxaGxMx).U75{^gcTf,uG^,)d{0)}\0) = q,T{q') + .(20) 

The leading contribution is determined by the following condensate, which 
is unknown: 

nq') = \mgcrx.Gx.yq). (21) 



■^The spectral density ImT{s) falls off quickly at large s in the chiral limit. This is an 
additional justification for the keeping only tt meson contribution. 



At the same time, both tt meson matrix elements which enter to this cor- 
relator are known - one of them: {Oldjp'j^igaxaG xau{x)\7T) is from (|T^) and 
another one: {Olu'j^^igaf^i^G ^u)d{0)}\n) is from (^) and can be found from 
PCAC Thus, we have finally the following very important equation: 



2 ^ . ^ „\ _ {qw(y\aGxaqf 



{qg (T^^G^yaxaGxaq) = pT^^ — *^^^^ 

Before we proceed with the numerical estimates, let me briefly formulate the 
result of the loop calculation for the correlator under consideration. After 
the "borelization" procedure the leading loop contribution to the T{M'^) 
is given by: T{M'^) = —-^M'^{qq), where M^ is Borel parameter which is 
determined by the next power corrections (which have the same small factor 
as/7r)and usually runs in the region M^ ~ 0.6 -r- IGeV"^. This term is 30 
times less than both contributions: ( PT| ) and vr meson residue (P^. This 
rate gives us some intuition about magnitude of the loop corrections to the 
eq.(P^. We neglect them in the rest of the paper. 



The factorization prescription for the condensate (|21|) leads to the follow- 
ing formula: 

K 3m^ 

{q9'(y,,G,,ax.Gx.q) = -j{9'G%G%){qq), K = ^^,^^ ^^ ^ ^ 3, (23) 

where we have introduced the coefficient of nonfactorizability K {K = 1 
if the factorization would work). The last equation follows from (^) and 
actually demonstrates that the factorization does not work. 

Few comments are in order. First of all, let us note that this result 
(the violation of factorization by a factor 3) is in a full agreement with the 
(absolutely independent) analysis [^ of the mixed vacuum condensates of 
the dimension seven. The result was based on the analysis of the heavy- 
light quark system and supports the present consideration. Besides that, 
let us note, that the VEV which appeared in the analysis []TT| was actually 



some combination of condensates {qg'^CF^yG^^axaGxaq) and (^S'^G^^G^j.g). 
We found for that combination more or less the same coefficient for the 
nonfactorizability i^ ~ 3. So, for the numerical estimates it is natural to 
assume that the factor i^ ~ 3 is the universal one for all condensates of this 
kind. 

As the second remark, we note that the phenomenon of non- factoriz- 
ability for the mixed condensates is not a big surprise, because we faced the 
analogous phenomenon early |jl2[ for the four-fermion condensates with an 



"exotic" Lorentz structure. 

As the last remark, and probably, the most important one, we want to em- 
phasize that the eq.(^) can be considered as the recurrent relation between 
{qigcfxaGxaQ) and {qg"^ cr ^^G ^^a x^G xal) with the dimensional coefficient rri^. 
We can repeat this procedure for arbitrary n with result: 

{q{ig<yx.Gx„Tq) = {mlriqq), (24) 



which we expect is to be correct one with a high accuracy ~ fts/vr 0. The 
formula (^) is the main result of this section. First of all it demonstrates 
a very simple relation between condensates and gives very important phe- 
nomenological information for modelekonstructors of the QCD vacuum. Sec- 
ondly, bearing in mind that these condensates are directly related to the 
moments (fcj"), the formula (|2^ ) gives nontrivial information about nonper- 
turbative pion wave function ilj{k\,x). 

5. Constraints on the nonperturbative wave function ^^(A;^,x). 

Now we are going to consider some applications of the obtained results. 
First of all let us consider the pseudoscalar wf and its moments. From eqs. 
(H,^ we have the following ratio for {k'i) p / {kj_)p: 

(^^^%!G^^^ (25) 

{kl)l 9mt ^ ^ 

where we assumed the universality for the nonfactorazibility factor i^ ^ 3, 
as explained above. Let us remark, that we know the condensate (|22|), which 
enters to the expression for {k\)p "almost" exactly. However, we do not 
possess such precise information for the second VEV which contributes to 
{k'\_)p with the same weight. In spite of this fact we expect a good accuracy 
for the ratio, because both contributions go with the same signs, and because 
some independent combination of these VEVs leads to the same coefficient 
i^ ^ 3, see above. From the pure theoretical point of view, by considering 
the limit d ^ oo, we can argue that there is not any cancellation between 
different terms (most dangerous thing which could happen!). 

This ratio actually is some expression of the fluctuations of the momen- 
tum. The result (^) demonstrates that the distribution function for the 
pseudoscalar function is rather compact in the momentum space in a big 
contrast with the analogous distribution for the leading twist wave function. 



For the axial wf from the formulae (|T8|,|T9|) we have the following ratio 
(instead of eq.(p5|)): 

which tells us that the fluctuations of the transverse momentum are much 
larger for the axial wf than for the pseudoscalar wf 0. 

■^The analogous result for the {q{G\a)"'q) looks as follows: {qiG^^C^^y- q) — {m\)^{qq), 
but with unknown coefficient m\. 

^Let us note, that the same ratio was examined in the recent preprint |13|| by the 

method of 0] with the result -W^T = ^ which numerically is not far away from our 
estimates (pq). However from our opinion, the corresponding sum rules are not suitable 
for the calculations of (fc")yi. We see at least few reasons for that. First of all, the 



It is interesting to note, that the analysis of the distribution function for 
the longitudinal momentum (^^") [Q shows the same phenomenon: deviation 
of longitudinal moments from their asymptotic values is large for the axial 
wf and very modest for the pseudoscalar wf. 

In principle we are ready to model ip{kj_,x). However, as is known, the 
knowledge of a finite number of moments is not sufficient to completely de- 
termine the ilj{kj_,x). The behavior of the asymptotically distant terms is 
a very important thing as well. Let me remind you of an example of the 
well-known case for the longitudinal distribution which determined by the 
0(^). If (^'") would behave as ~ 1/m at large m, it would mean that the 
integrand J_i <i^^™'0(^) ~ 1/m implying that 0(^ — > ±1) — ^ constant. At the 
same time, the very likely assumption that the vr-meson fills a finite duality 
interval in the dispersion relation at any m means that J_i c?^^™0(^) ~ l/rn^ 
0. It unambiguously implies the behavior 

.1 0(^^±l)-.(l-^2^. 

A simple model wave function which meets this requirement as well as pos- 
sesses a large value of the lowest moment, (^^) ~ 0.4 has been proposed in 
ref. i: 

1 "i 

0(o = -(i-ae', (0 = i, (a = 0.43 (27) 

We want to emphasize that the "endpoint" behavior (which corresponds to 
asymptotically distant terms) was crucial in this analysis. 

In terms of iplk"^,^ = 2x — 1) this constraints can be rewritten in the 
following form: J dk'j_ilj{k'j_, C, -^ ±1) ~ (1 — ^^). The analogous assumption, 
that the vr-meson fills a finite duality interval in the corresponding dispersion 
relation at any n, where n is related to n- transverse moment (k"^) gives the 
following constraint: 



|t/^i^ix^i,e-±i)~(i-e 



This constraint is extremely important and implies that the k^ dependence of 
the 'il){k']_.,^ = 2x — 1) comes exclusively in the combination A;^/(l — ,^^) at 
^ —>■ ±1. The byproduct of this constraint can be formulated as follows. The 
standard assumption on factorizability of the ip{^±yO = ^"(^1)0(0 does 
contradict to the very general property of the theory formulated above. 

corresponding spectral density increases very quickly ~ s" with s and there is no reason 
for IT meson saturation. As a consequence of it, there is a very strong dependence on 
So which is far away {So ~ OAGeV^) from the standard value 5*0 — IGeV^. By some 
kinematical reasons, the 4-quark condensate, which used to be the most important one, 
does not contribute. It clearly means, that the next power corrections, proportional to 
{q{igcr\aGxc)'^q){qq) become very important. Thus, our point is that the corresponding 
sum rules do not work and moments (fc^") can not be extracted in such way. 



The next constraints come from the calculation {kY)Ai (|T8|) and (fc^), 

• Sjdkl j\ d^kl, = 1, ikl) = ^ ^ O.lGeV^ {ki) » {kiy. 

These constraints actually give the general scale of the ip{k'j_, C,). Besides that, 
the large fluctuations of the momentum mean, first of all, that the distribu- 
tion of the transverse quark momentum is wide and it is not concentrated 
about kj_ ~ 0. 

Our last constraint comes from the analysis of the asymptotically distant 
terms. We do not know the behavior of the asymptotically large moments 
(fc^") ai n —>■ oo exactly. But we do know the dependence on n in the limit 
when the space-time dimension rf — > oo, (|^)Q. We expect (and this is the 
main assumption) that in the real world the functional dependence will not 
be changed. Thus, for the asymptotically large n we expect the following 
behavior: 

where at the last stage we have used the eq.(^), which is suppose to be 
valid for the arbitrary n. It is very important that the right hand side of this 
equation is constant and not zero. We are not going to use this constant in 
our model as input parameter; the only fact we need, that this is not zero. 
Let us note, that any mild function on n, like 1/n or even exp(n) on the right 
hand side can not be ruled out. Theonly effect it brings, is somejescale of 
the dimensional parameter: m^ — >■ ttiq = ^^rn^ = rn^ or rn^ — > ttiq = emg, 
which anyhow, is unknown. The same constraint can be obtained from the 
dispersion relation with the assumption that the tt meson fills a finite duality 
interval 5*0 at any n, like in the analysis (•!). In this case, instead of (•4) 
we get 



fl{kl^) 



35o"+'(2n + 2)!! 8"^+^ 
87r2(2n + 3)!! n 



which can be reduced to the previous one. 

Very important consequence of this constraint can be easily seen if we 
rewrite it in the following form: 

•4 j dkl^l^i^^Mlr ^ I. 

which means that the 'ip{jzFi) ^ ^( w2 — S) for sufficiently large k^. We 
have to pause here in order to explain the general idea of the Wilson operator 

^We derived this formula for pseudoscalar function only, but absolutely the same for- 
mula takes place for the axial wf as well. The derivation is the same. 
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expansion (OPE) in the given context and the term "sufficiently large fc^" 
in particular. 

As is known, all VEVs within OPE are defined in such a way, that all 
gluon's and quark's virtualities, smaller than some parameter /x (point of 
normalization) are hidden in the definition of the "nonperturbative vacuum 
matrix elements". All virtualities larger than that should be taken into ac- 
count perturbatively, see for detail discussion of this problem. At the 
same time, from the exact PCAC relations, like (§) or(|^) it is clear that mo- 
ments of the wave functions are related to condensates, defined as explained 
above. Thus, all transverse moments are defined in the same way as conden- 
sates do. This is, actually, the definition of tlie nonperturbative wave 
function ijj{k']_,^), through its moments which can be expressed in terms of 
nonperturbative vacuum condensates. 

Now it clear, what we mean by term "sufficiently large kj_". By that 
we mean the largest virtuality which have been taken into account in con- 
struction oi iplk"^,^), or (what is the same) in the definition of condensates. 
The corresponding dimensional parameter can be expressed in terms of some 
condensate (see bellow) and it depends on fi only logarithmically, like con- 
densates do. The important consequence of the definition can be formulated 
in the following way: The nonperturbative wf even for sufficiently large k'j_ 
does not behave like l/k'j_ as is frequently assumed. In order to understand 
this statement let us imagine that we calculate the gluon condensate (G^^). 
It is clear that we will not substitute the gluon propagator in the form 1/fc^ in 
this calculation, because, by definition, the perturbative contribution should 
be subtracted in the calculation of the nonperturbative condensate. Have in 
mind that our definition of the nonperturbative wf is formulated in terms 
of the vacuum condensates, it is clear that the same statement is concerned 
to the nonperturbative wave function as well and its behavior has nothing to 
do with l/kj_. 

We want to model the simplest version of the wave function which meets 
all requirements formulated above. First of all, as was explained, the kj_ 

dependence comes only in the combination ipi^iPi) in order to meet the re- 
quirement (•2). In order to satisfy the constraint (•4), we have to assume 

that at sufficiently large k^ we have a S{yzci — S) -like function, whose mo- 
ments correctly reproduce the constant on the right hand side of the equation; 
the S is some input dimensional parameter, which will be expressed in terms 
of (kl). 

Our next step is to satisfy the constraint (•S). It is clear that the 6 
function proposed above does not provide a noticeable fluctuations of the 
momentum In order to meet this requirement, we have to spread out the 
distribution function between fc^^ ~ and k'j_ ^ S in such a way, that the 
overall area will be the same, but moments should satisfy to this requirement. 
In principle, it can be done in arbitrary way. The simplest way to make the 
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wf wider is to put another 5 function at kj_ = 0. With these remarks in 
mind we propose the following "two- hump" {again f) nonperturbative wave 
function which meets all requirements discussed above: 

^(kl, = {A5i^, -S) + B6i^)Me - ^) + ^] 

jdklf^diij{kl,i) = l, ^ + B = ^ (28) 

m ^jdki^{ki,o = f (1 - erne - ^) + J] 

(fci)^5(fcl)2^A = 7/8, g = l, {kl) = ^-^^S=^-^{kl)c^0.8GeV' 

6b 2 

Few comments are in order. We put the common factor [^'(^^ — |) + |] in 
the front of the formula in order to reproduce the light cone 0(^) function 
with arbitrary (^^). For 51 = it corresponds to the asymptotic wf : 0(^) = 
3/4(1 - ^2). For ^ = 1 we reproduce the (j){Ocz = 15/4(1 - ^2)^2 ^-^^-^ 
(^^) ^ 0.43. We still believe in large value for the moment (^^) in spite of 
the criticism from the ref. |14|, who found much less value for (^^). This is 
not the place to discuss this question in more detail, but I want to make a 
comment that there are few questions to be answered ( theoretical, as well as 



phenomenological ones) before the approach advocated by authors of ref. [^ 
can be considered as is well defined and based on solid ground. At the same 
time, the people who prefer to use a smaller value for (^^), can make their 



own choice by changing the parameter g in the eq. (|28|) . Let me remind the 



relation (^^) = 1/5 + 8g/35 between (^^) and parameter g from the formula 
(P^) for doing so. It will not spoil the constraints discussed above. 

Let me emphasize, that the formula (^) is not destined for the precise 
fitting of the experimental data and it is given as illustration only. By the 
physical reasons it is clear that the 6 functions should be spread out in 
some way. This procedure is definitely model dependent. We only want to 
attract attention on the asymmetric form of the wf in the k± space. This 
qualitative property of the wf comes from the large magnitude of the ratio 
r = {ki)/{kir. 

As the last remark, let me point out, that the description of exclusive 
reactions at experimentally accessible momentum transfers {IGeV"^ < Q^ < 
306*6^^) and the problem of extracting some information about nonpertur- 
bative wf are two different problems. It is was advocated by Radyushkin 
and collaborators [1^, [|1^ and Isgur, Llewellyn Smith [|1^ for a long time 
that the asymptotically leading contributions do not describe the experimen- 



tal data at moderate Q^, for recent review of this subject, see [18|. Some 
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confirmation of this statement came recently from |^, where it was shown 
that the inclusion of the intrinsic k'j_- dependence and Sudakov suppression 
leads to the self-consistent calculation, but the obtained magnitude is too 
small with respect to the data even if (pcz asymmetric function is used. It 
is very likely, that the "soft" contributions play an important role in this 
region. 

It turns out, that the very unusual shape of the ip{kj_, ^) described above 
supports this idea and can imitate the behavior of the asymptotically leading 
contribution at the intermediate momentum transfer. Such a mechanism, if 
it is correct, would be an explanation of the phenomenological success of the 
dimensional counting rules at a very modest Q^. We are going to discuss this 
question in more detail somewhere else [|19] . 

6. Conclusion. 

Let me formulate the main results of this paper. 

First of all, we used the standard definition of the nonperturbative wf 
through its moments and we expressed these moments in terms of the vacuum 
condensates. From this construction it becomes clear that the definition of 
the vacuum condensates (within OPE), with subtraction of the perturbative 
contributions, and the definition of the nonperturbative wf is one and the 
same problem. 

Secondly, we formulated some constraints on the nonperturbative wf; 
they have very general origin, and thus, they should be satisfied for any 
model. There is no reason to repeat these constraints (•! — •4) from the pre- 
vious section again; the only comment we want to make here is the following: 
the standard assumption on factorizability of the wf : ip{kj_, x) = ilj{k\)(j){x) 
can not be matched with these constraints. 

The last result which deserves to be mentioned here and which probably 
is interesting by its own, regardless to the question of nonperturbative wf, 
looks as follows {q{iga\aGxa-Yq) = {jnlY{qq). I am not aware of any model 
of the QCD vacuum, which satisfies to this "almost" exact relation between 
mixed vacuum condensates. 
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